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Abstract 

The electroproduction of J/ip and ip(2S) mesons is studied in elastic, quasi-elastic and 
inclusive reactions for four momentum transfers 2 < Q 2 < 80 GeV 2 and photon-proton 
centre of mass energies 25 < W < 180 GeV. The data were taken with the HI detector 
at the electron proton collider HERA in the years 1995 to 1997. The total virtual photon- 
proton cross section for elastic J/ip production is measured as a function of Q 2 and W. 
The dependence of the production rates on the square of the momentum transfer from the 
proton (t) is extracted. Decay angular distributions are analysed and the ratio of the longi- 
tudinal and transverse cross sections is derived. The ratio of the cross sections for quasi- 
elastic ip(2S) and J/ip meson production is measured as a function of Q 2 . The results are 
discussed in terms of theoretical models based upon perturbative QCD. Differential cross 
sections for inclusive and inelastic production of J/ip mesons are determined and predic- 
tions within two theoretical frameworks are compared with the data, the non-relativistic 
QCD factorization approach including colour octet and colour singlet contributions, and 
the model of Soft Colour Interactions. 
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1 Introduction 



The high energy electron-proton collider HERA has renewed the interest in the study of light 
and heavy vector mesons produced in processes with quasi real and virtual photon exchange. 
Several production mechanisms valid in limited kinematic regions have been discussed in the 
literature for such processes and a unified picture is not available. The topics of the present paper 
are studies of elastic and inelastic production of J/ip mesons and of quasi-elastic production of 
ip(2S) mesons for four momentum transfers 2 < Q 2 < 80 GeV 2 and photon-proton centre of 
mass energies 25 < W < 180 GeV. 

Leptoproduction of J/ip mesons has previously been studied in several fixed target experiments 
and at HERA in different kinematic regions [1-3]. In photoproduction, which corresponds to the 
limit Q 2 ~ 0, and at low and medium Q 2 the production of J/ip mesons, e + p — > e + J/ip + X, 
is found to be dominated by processes where the hadronic system X is either a proton (elastic 
J/ip production) or has a low mass M X - These processes show characteristics of diffraction 
as observed in hadron-hadron interactions and in photoproduction of light vector mesons at 
lower energies. However, the experiments HI and ZEUS have found that at HERA energies the 
dependence of the elastic J/tp cross section on W in the photoproduction limit is steeper than 
measured in soft diffractive processes B|j. 

In recent years it has been demonstrated that elastic photo- and electroproduction of J/ip 
mesons (Fig. |T]a, c) can be calculated within perturbative QCD (pQCD) [6-9]. In these cal- 
culations the elastic cross section is related to the square of the gluon density in the proton 
and the fast rise of elastic J/ip production with W reflects the increase of the gluon density 
in the proton at low values of Bjorken x JT0|] . According to these models elastic J/ip meson 
production consequently offers a sensitive way to probe the gluon density. Further predictions 
of such models concern, for example, the fraction of longitudinally polarized J / ip mesons, the 
dependences of the slope of the t distribution (t is the square of the momentum transfer from 
the proton) and of the ratio of ip(2S) to J/ip meson production on kinematic variables. The 
latter ratio is also predicted in an approach based upon colour dipole phenomenology [ |TT| ] . 

Inelastic J/ip production, which can be described by the formation of cc states via boson gluon 
fusion (Fig. [j]d, e), was previously studied in fixed target experiments [jl], |12[] and was analysed 
in the framework of the Colour Singlet Model [|T3|] . At HERA inelastic J/ip production has 
been analysed in photoproduction [Q, [14|] and was successfully described by predictions of the 
Colour Singlet Model in next-to-leading order JT3|]. On the other hand, measurements of the 
production of J/ip mesons in hadronic collisions [ [T6| ] have shown that the Colour Singlet Model 
cannot account for the observed cross section. A good description of these data can, however, 
be achieved using a factorization approach in the framework of non-relativistic QCD (NRQCD) 
[jT7p, where also colour octet states contribute. These colour octet processes are also expected 
to contribute in electroproduction. First analyses of HERA data showed that the color octet 
contribution in the photoproduction regime is less than expected [Jl8j], In order to shed further 
light on the production process we present a fully inclusive analysis of J / ip meson production in 
the range 0.2 < z < 1, where z is the ratio of the energies of the J/ip and the exchanged photon 
in the proton rest frame. In addition we extract inelastic cross sections in the same z range. 
The data are compared to calculations of lepton proton scattering in leading order performed 
in the NRQCD formalism [ |19T ] and to a phenomenological model incorporating Soft Colour 
Interactions [|2(]] in the Monte Carlo generator AROMA ^TJ]. 
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Figure 1: Charmonium production mechanisms: a) elastic and b) proton dissociative production 
via pomeron exchange; c) elastic production via two-gluon exchange; d) leading order diagram 
in the Colour Singlet Model (the cc pair is produced in a colour singlet state); e) leading order 
Colour Octet Model (the cc pair is produced in a colour octet state). 

The paper is organised as follows: after a discussion of different charmonium production models 
relevant for the present analysis and a description of the event selection, the total and differential 
cross sections for the elastic reaction e + p — > e + J/ip + p are presented with an extended 
kinematic reach compared to our previous measurement [Q] and with statistics increased by 
an order of magnitude. We then report on the first measurement of ip(2S) production in deep 
inelastic scattering at HERA and extract the ratio of the ip(2S) to the J/ip meson production 
cross section as a function of Q 2 . Finally differential and total cross sections for inclusive and 
inelastic J/ip production are presented. 



2 Models and Phenomenology 

The experimental distinction between the various J/ip production mechanisms discussed in the 
literature is not unambiguous and the following terminology will be adopted here. The process 

e + p -> e + J/ip + X 

will be called "elastic" if X is a proton. Since the proton is in general not observed we use the 
term "quasi-elastic" for events in which only the tracks of the J/ip decay leptons are present 
in the main detector. This data sample comprises in addition to elastic events those in which 
the proton is diffractively excited into a system X dominantly of low mass, which subsequently 
dissociates (Fig. [l]b). The decay or fragmentation products of this low mass system in general 
escape detection in the main detector. If a high mass system X is produced the emerging 
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hadrons are usually detected and the process is called "inelastic". The term "inclusive" is used 
if only the presence of a J/tp is demanded irrespective of the production and detection of other 
particles. 

Elastic Charmonium Production Elastic photoproduction of light vector mesons, p, uj and 
0, is characterized by a weak dependence of the cross section on the photon-proton centre of 
mass energy W and by a diffractive peak, i.e. small scattering angle of the vector meson with 
respect to the incident photon direction. This behaviour is well described by vector meson 
dominance and Regge theory in terms of soft pomeron exchange (Fig. |TJa). However, the cross 
section for elastic production of J/ip mesons by quasi-real photons (Q 2 ~ 0) at HERA is 
observed to rise steeply with W. Parameterizing the dependence as W 5 , S is measured to 
be of order 1 for J/ip mesons [Q, while light vector mesons show an energy dependence 
compatible with expectations from pomeron exchange in soft hadronic processes corresponding 
to 5 ~ 0.22 - 0.32 (see ©). 

Modifications of the simple soft pomeron exchange model were subsequently proposed to de- 
scribe the HERA data [ |2"3| ] . Alternatively an approach in the framework of pQCD [7], |9|] 
was pursued. In the models based on pQCD the interaction between the proton and the cc pair 
is mediated by a system of two gluons (Fig. [I]c) or a gluon ladder and the fast increase of the 
cross section is related to the rise of the gluon density in the proton at small values of Bjorken 
x. Since the gluon density enters the cross section quadratically the sensitivity is large. In these 
calculations the scattering amplitude is obtained from the convolution of three contributions 
which are characterized by different time scales: the fluctuation of the (virtual) photon into a 
cc pair, the scattering of this hadronic system on the proton, and the formation of the final state 
vector meson. 

In contrast to photoproduction where the charm quark mass offers the only hard scale (at low 
values of \t\), electroproduction of <//?/> mesons has an additional scale, Q 2 , and at high Q 2 the 
predictions of perturbation theory are expected to become more reliable. Electroproduction of 
heavy vector mesons within pQCD was recently studied in great detail by Frankfurt et al. [|9|]. 
Important corrections were found concerning, for example, the choice of the scale at which the 
gluon density is probed, concerning the gluon distribution and the wave function for the vector 
meson and the importance of corrections due to Fermi motion of the quarks within the vector 
meson. 

Inelastic Charmonium Production In the Colour Singlet Model photo- or electroproduction 
of J/ij) mesons is assumed to proceed via boson gluon fusion into cc pairs which emerge in a 
colour singlet state due to the emission of an additional hard gluon (Fig. [T|d). The failure of the 
Colour Singlet Model in describing hadroproduction of quarkonia at large transverse momenta 
led to new approaches for the description of J/tp production which include contributions from 
cc pairs in colour octet states. The approach by Bodwin, Braaten and Lepage (BBL) [ |T7| ] based 
on a factorization approach in NRQCD was suggested to describe the large hadroproduction 
rates for J/ip and ip(2S) production at large p t measured at the Tevatron p% 
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When applied to J/ip electroproduction [|19"|], the cross section in the BBL formalism can be 
expressed as: 



a(e + p^e + J/tp + X) = c n(e +P e + cc[n) + X)(0^), (1) 

rt 

where cc[n] denotes an intermediate cc pair in a definite colour, spin and angular momentum 
state n. For each n, the cross section factorizes into a short distance part c n calculable in a 
perturbative QCD expansion in the strong coupling parameter a s and a long distance matrix 
element (On^) representing the probability for the cc[n) pair to evolve into a colour singlet 
J/tp meson and additional soft gluons. The long distance matrix elements are not calculable 
in perturbation theory and have to be determined experimentally or by lattice calculations, but 
they are thought to be universal. 

The relative importance of the terms in equation (Jl]) is determined by NRQCD scaling laws with 
respect to v, the typical relative velocity of the charm quarks in the cc system. In contrast to the 
Colour Singlet Model in which all c n not corresponding to colour singlet cc states are neglected, 
the BBL formalism includes states where the cc system is in a colour octet state (Fig. |I]e); 
therefore it is often called the "Colour Octet Model". However, colour octet contributions are 
suppressed by powers of v 2 . Since v is small, (v 2 ) ~ 0.3 []T7[], they only become important 
when the corresponding short distance coefficients c n are large. In the limit v — > the Colour 
Singlet Model is restored. 



Soft Colour Interactions The model of Soft Colour Interactions (SCI) was originally devel- 
oped as an alternative to Regge phenomenology and pomeron exchange to describe diffractive 
scattering at HERA [fZ5p. The model was successfully applied to quarkonia production at the 
Tevatron 



For electroproduction the model was implemented in the Monte Carlo program AROMA P1Q 
which generates cc pairs via photon gluon fusion according to the leading order matrix ele- 
ments approximating higher orders by parton showers. At a scale below the cut-off for pQCD 
additional interactions take place: quarks and gluons generated in the hard process interact non- 
perturbatively with the partons of the proton remnant; the latter are also allowed to interact with 
each other. In these soft colour interactions the momenta of the partons are not affected, only 
their colour states may change, which leads to a modification of the hadronic final state. 

In such a model the conversion of a primary cc pair - which is in a colour octet state - into an 
observable colour singlet state such as J/ip, Xc, etc. may occur if the mass corresponds to the 
mass of the produced charmonium particle. The probability for this to happen is not constant as, 
for example, assumed in the Colour Evaporation Model [ ^6| ] but depends on the final partonic 
state and, in the Lund string model, on the string configuration. 



3 Detector, Event Selection, Kinematics and Simulations 

The data presented here correspond to an integrated luminosity of 27.3 ± 0.4 pb _1 . They were 
collected in the years 1995 to 1997 using the HI detector which is described in detail in p7j]. 
HERA was operated with 27.5 GeV positrons and 820 GeV protons. 
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3.1 Detector and Event Selection 



J /if) mesons are detected via the decays J /if) — > /i + \i and J /if) — > e + e with branching 
fractions of 6.01 (6.02) ± 0.19%, respectively g§]. For if)(2S) mesons the decay ip(2S) -> 
J/ipn + n~ is used (branching ratio 30.2 ± 1.9% [f£BJ]) with the subsequent decay of the J/-0 
into or e + e _ . The criteria for the data selection are summarised in Table |T]; further details 
of the analysis can be found in [E9p. 



In the Q 2 range studied here, the scattered positron is identified by its energy deposition in the 
backward electromagnetic calorimeter SpaCal pO] ] situated 152 cm backward from the nominal 
interaction point of the beamsQ. The SpaCal covers the polar angles 155° < 9 < 178° and has an 
energy resolution of cr E /E ~ 7.5%/ v /£/GeV©2.5%. Aminimal energy deposition of 14 GeV 
is required and cuts are applied to the cluster position and cluster shape in order to ensure high 
trigger efficiency and a good quality positron measurement. To keep acceptance corrections 
small, the Q 2 range is limited to 2 < Q 2 < 80 GeV 2 . A drift chamber (BDC) in front of the 
SpaCal is used to reconstruct the polar angle 9 e of the scattered positron in combination with 
the interaction vertex. 

The decay leptons of the J/if> meson are detected in the central tracking detector (CTD), con- 
sisting mainly of two coaxial cylindrical drift chambers, which have a length of 2.2 m and outer 
radii of 0.45 m and 0.84 m. The charged particle momentum component transverse to the beam 
direction is measured in these chambers by the track curvature in the 1.15 T magnetic field gen- 
erated by the superconducting solenoid, with the field lines directed along the beam axis. Two 
polygonal drift chambers with wires perpendicular to the beam direction, which are located re- 
spectively at the inner radii of the two chambers, are used to improve the measurement of the 
particle polar angle. The tracking system is complemented in the forward direction by a set of 
drift chambers with wires perpendicular to the beam direction which allow particle detection 
for polar angles 9^7°. Multiwire proportional chambers serve for triggering purposes. 

In the present analyses, two oppositely charged tracks with transverse momenta p t larger than 
0.1 GeV are required to be reconstructed in the CTD with polar angles in the range 20° < 9 < 
160° where the detection efficiency is high. For each event, the vertex position in z is determined 
using tracks reconstructed in the CTD. To suppress background from interactions of the beam 
with residual gas in the beam pipe, the vertex must be reconstructed within 40 cm from the 
nominal interaction point corresponding to 3.5 times the width of the vertex distribution. 

For the study of elastic J /if) production to suppress background from inelastic reactions no 
track, except a possible track from the scattered positron, is allowed to be present in addition to 
the tracks from the two decay leptons. In the if>(2S) analysis exactly two tracks, assumed to be 
pions, with transverse momenta above 0.12 GeV and opposite charge are required in addition to 
the two decay leptons (an additional track from the scattered positron is allowed). In the study 
of inclusive J/if> production no requirement on the track multiplicity is imposed. An inelastic 
data set is defined by the additional requirement of an energy deposition of Ef w d > 5 GeV in 
the forward region of the liquid argon (LAr) calorimeter at polar angles 9 < 20°. 

The J /if) decay leptons are identified by the LAr calorimeter surrounding the tracking detec- 
tors and situated inside the solenoid. The LAr calorimeter is segmented into electromagnetic 

'HI uses a right handed coordinate system, the forward {+z) direction, with respect to which polar angles are 
measured, is defined as that of the incident proton beam, the backward direction (— z) is that of the positron beam. 
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I. Quasi-Elastic J/-0 CTD-CTD 40 < W < 160 GeV, 2 < Q 2 < 80 GeV 2 



Tracks 


Reconstructed event vertex with \z vtx — z nom \ < 40 cm 
Exactly 2 tracks in CTD 2 

Opposite charges, 20° < 9 < 160°, p t > 0.1 GeV 


Decay leptons 


> 1 /i identified in LAr Cal. (LAr) or Central Muon Detector (CMD) or 
2 e identified in LAr Calorimeter 


Other 


Forward untagged: Ej°2 r < 1 GeV and N PRT = and N FMD < 2 
Forward tagged: Ef° Ar > 1 GeV or N PRT > or N FM d > 2 


II. Quasi-Elastic J/-0 FMD-CTD (FMD-FMD) 


25 < W < 40 GeV, 2 < Q 2 < 6 GeV 2 


Decay fi 


1 n in FMD and 1 // in CTD+LAr/CMD as in I. or 

2 fx in FMD 
Opposite charges 


Other 


No tracks except those associated with the decay muons 2 


III. Quasi-elastic ^(2S) 40 < W < 180 GeV, 1 < Q 2 < 80 GeV 2 


Decay particles 


Reconstructed event vertex with \z vtx — z nom \ < 40 cm 
Exactly 4 tracks in CTD 2 

2 /i identified in LAr or CMD or 2 e identified in LAr 

Opposite charges of the two leptons and of the two additional tracks (pions) 

20° < 9 < 160°, p t > 800 MeV (leptons), p t > 120 MeV (tt+, tt") 


IV. Inclusive and Inelastic J/ip 40 < W < 180 GeV 


, 2 < Q 2 < 80 GeV 2 


Tracks 


Reconstructed event vertex with 




< 40 cm 


Decay leptons 


2 /i identified in LAr or CMD or 2 e identified in LAr 
Opposite charges, 20° < 9 < 160° a.ndp t > 800 MeV 


Inelastic selection 


E fwd > 5 GeV 


General 


Scattered e + 


E e > 14 GeV identified in SpaCal 


Final state 


E(E-Pz) >45 GeV 



2 Any additional track associated with the scattered e+ is not considered here. 



Table 1: Summary of selection criteria for the different data sets. 
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and hadronic sections and covers the polar angular range 4° < 6 < 154° with full azimuthal 
coverage. Muons which are identified as minimum ionizing particles in the LAr calorimeter 
can in addition be identified by track segments reconstructed in the instrumented iron return 
yoke (central muon detector CMD, 4° < 9 < 171°) and in the forward muon detector (FMD, 
3° < 9 < 17°). The FMD provides track segments in front of and behind a toroidal magnet 
with a field of B = 1.5 — 1.75 T, thus allowing a determination of the muon momentum with a 
precision of about 25%. 

The lepton selection criteria vary for the different data sets, depending on the amount of back- 
ground. For the quasi-elastic event selection, the identification of one muon or two electrons 
is required. For the ip(2S) and the inclusive data sets two identified leptons with transverse 
momenta p t > 0.8 GeV are required. 

A dedicated analysis has been carried out to extend the analysis of elastic J/ip production 
towards small W (25 < W < 40 GeV) using the FMD ("Low W analysis"). Two data sets 
are selected. The events of the first set which are required to have one decay muon in the 
FMD and the other one in the CTD (FMD-CTD) are used for a cross section determination. 
In the second data set both muons are measured in the FMD (FMD-FMD) with an additional 
loose vertex requirement to suppress muons originating from the proton beam halo. The latter 
set (20 < W < 35 GeV) serves as a control sample for the FMD efficiency which was also 
determined from a larger photoproduction J/ip sample. Due to statistical limitations, the low 
W analysis is restricted to the low Q 2 (2 < Q 2 < 6 GeV 2 ) region. 

The triggers in all analyses require a total energy deposition in the SpaCal above a threshold 
of 6 — 12 GeV. The value of the threshold depends on the topology of the energy deposition 
and on the presence of additional requirements, such as signals from the central drift chambers 
and/or the multiwire proportional chambers. 

In order to minimise the effects of QED radiation in the initial state, the difference between 
the total energy and the total longitudinal momentum ^2(E — p z ) reconstructed in the event is 
required to be larger than 45 GeV. If no particle, in particular a radiated photon, has escaped 
detection in the backward direction, the value of ^2(E — p z ) is twice the incident positron 
energy, i.e. 55 GeV. 

Forward Region After requiring exactly two tracks corresponding to the J/ip decay leptons 
the event sample includes two main contributions: elastic events and events with proton dissoci- 
ation. It is possible to identify most of the proton dissociation events with the components of the 
detector in the forward region, namely the forward part of the LAr calorimeter, 4° < 9 < 10°, 
the forward muon detector FMD and the proton remnant tagger (PRT, an array of scintillators 
24 m downstream of the interaction point, 0.06° < 9 < 0.17°). When particles from the diffrac- 
tively excited system interact with the material in the beam pipe or with the collimators, the 
interaction products can be detected in these forward detectors. Events are "tagged" as candi- 
dates for proton dissociation by the presence of a cluster with energy E LAr larger than 1 GeV 
at an angle 9 < 10° in the LAr calorimeter, or by at least 2 pairs of hits in the first three layers 
of the FMD (Nfmd > 2), or by at least one hit in the proton remnant tagger (Nprt > 0). The 
forward detectors are sensitive to M x £1.6 GeV. 
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3.2 Kinematics 



The kinematics for charmonium production is described with the standard variables used for 
deep inelastic interactions, namely the square of the ep centre of mass energy, s = (p + k) 2 , 
Q 2 = —q 2 and W = (p+q) 2 , where k, p and q are the four- momenta of the incident positron and 
proton and of the virtual photon. In addition, the scaled energy transfer y = p ■ q/p ■ k (energy 
fraction transferred from the positron to the hadronic final state in the proton rest frame) and 
the Bjorken variable x = Q 2 /2p ■ q are used. Neglecting the positron and proton masses the 
following relations hold: Q 2 = xys and W 2 = ys — Q 2 . 

In the case of elastic J/ip and quasi-elastic ip(2S) production the kinematic variables are re- 
constructed with the "double angle" method [|31"|], where Q 2 and y are computed using the polar 
angles 9 and 7 of the positron and of the vector meson in the HERA laboratory frame of refer- 
ence, which are well measured: 



_, 2 sin 7(1 + cos 1 
sin 7 + sin 9 — sin(7 + 9) 



Q 2 = 4£g - • V , (2) 



_ sin 0(1 - cos 7) ^ 
^ sin 7 + sin 6* — sin(7 + 9) ' 

E denotes the energy of the incident positron. The momentum components of the J /if) and the 
ip(2S) mesons are obtained from their measured decay products. 

Since the fractional energy loss of the proton is negligible, the absolute value of the four mo- 
mentum transfer t is given to a good approximation by the following relation^: 

1*1 - (Ptp? = (Pte+PtvY, (4) 

where p tp , pt e and p tv are, respectively, the momentum components transverse to the beam 
direction of the final state proton, positron^ and vector meson. The resolution obtained from the 
Monte Carlo simulation for the reconstruction of W ranges from 4 to 9% depending on Q 2 . For 
Q 2 it is about 2%, and for t on average 0.10 GeV 2 . 

The variable (E ~ Vz) is computed as: 

J2( E - Vz) = {E e + E v ) - {p ze + p zv ), (5) 

E e and E v being the measured energies of the scattered positron and of the vector meson, and 
p ze and p zv their momentum components parallel to the beam direction. 

For inclusive J/ r tp production, the "eS" method [ p2| ] is used to reconstruct the event kinematics, 
which combines the measurement of both the scattered positron and the full hadronic final state 
to obtain good resolution over the entire kinematic region. The variable Q 2 is reconstructed 

3 The lowest \t\ value kinematically allowed, t m in — (Q 2 + m 2 /) 2 rrip /y 2 s 2 , is negligibly small. 
4 The momentum of the scattered positron is here computed from Q 2 and y, which provides better precision 
than the direct SpaCal energy measurement. 
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from the scattered positron. For the calculation of y and the elasticity z = (p^ •p)/(q-p), where 
PtJi denotes the J/ip four-momentum, the observed final state is used in addition. Thus 

y= E *7f and ,= < £ -*■>* , (6) 

where the sums run over all particles observed in the final state, but excluding the scattered 
positron in those indicated by "had". For the calculation of the sums in equation @ a combi- 
nation of tracks reconstructed in the CTD and cells in the LAr and SpaCal calorimeters is used. 
The resolution is good (2 — 5%) for the variables Q 2 , p\^ and y*, the rapidity of the J/ip in the 
Yp centre of mass frame. For z (W) the resolution is 2% (3%) for z > 0.9 and on average 17% 
(9%) for z < 0.9. 



3.3 Monte Carlo Simulations 

To take account of detector acceptance and efficiencies, smearing effects, losses due to the 
selection criteria, and remaining backgrounds, corrections are applied to the data using Monte 
Carlo simulations. The HI detector response is simulated in detail, and the simulated events are 
passed through the same reconstruction and analysis chain as the data. 

The correct description of the data by the simulation has been checked extensively by inde- 
pendent measurements and was adjusted where necessary. In particular, the trigger efficiencies 
have been determined using independent data sets, the efficiency of the central drift chambers 
has been measured using cosmic ray muons, and the lepton identification probabilities have 
been determined using control samples in which the J/ip is reconstructed identifying only one 
or no lepton. Remaining differences between data and simulation are used to estimate system- 
atic uncertainties. 

Typical efficiencies are: lepton identification ~ 80% per lepton, track reconstruction ~ 96% per 
track, and identification of the scattered positron 99%. The total trigger efficiency is determined 
to be 97% on average. 

The following event generators are used: 



The DIFFVM program [ |33| ] is based on the Vector Meson Dominance Model and per- 
mits variation of the Q 2 , W and t dependences, as well as a variation of the value of 
R = gl/ot- In addition to the elastic process, vector meson production with proton dis- 
sociation is simulated where the dependence of the cross section on the mass Mx of the 
dissociated hadronic state X is parameterized as 1/Mjp 6 '. High mass states are assumed 
to decay according to the Lund string model [[$4p. In the resonance domain, the mass 
distribution is modelled using measurements from target dissociation on deuterium p5|], 
and resonance decays are described using their known branching ratios. 



EPJPSI [TJ^] implements inelastic J/ip production according to the Colour Singlet Model 
taking into account relativistic corrections and parton showers. 



The LPAIR generator [07p simulates QED electron- and muon-pair production, 77 



e + e and 77 — > , where the photons originate from the positron and proton respec- 
tively. Elastic and inelastic processes are simulated. 
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Radiative Corrections The measured cross sections are given in the QED Born approxi- 
mation for electron interactions. The effects of higher order processes - mainly initial state 
radiation - are estimated using the HECTOR program [|3"%|]. 

Radiative corrections for the double angle reconstruction (elastic J/ip and quasi-elastic ip(2S) 
production) are about 2— 3%, and are weakly dependent on Q 2 and W. A systematic uncertainty 
of 3% is obtained by variation of the Q 2 and W dependences of the 7*p cross section within 
the uncertainties of the measurement. For the reconstruction of kinematics according to the eS 
method (used in the inclusive J ftp analysis) the radiative corrections amount to 6 — 8%, and 
again are only weakly dependent on Q 2 and W. 



4 Elastic J I ^ Production 

The distributions of the invariant mass rriu for the selected events in the J /if) region with two 
tracks in the central region are presented in Fig. ||a and b, for the ^ and the e + e~ decay 
channels, respectively (data set CTD-CTD, I. in Table []]) . A clear signal is observed at 3.090 ± 
0.005 GeV, compatible with the nominal J/ip mass of 3.097 GeV [fZB]]. The peak width is 
compatible with the expectation obtained from the detector simulation. The mass spectra for 
the events in the low W analysis where one or both muons are reconstructed in the forward 
muon detector (FMD) are shown in Fig. |a and b (data sets FMD-FMD and FMD-CTD, II. in 
Table [I]). 

For the events with two tracks in the central detector the J/ip signal region is defined by the 
condition \mu — < 250 MeV, where is the nominal J/ip mass. The non-resonant 
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Figure 2: Mass spectra for events of the quasi-elastic J ftp selection: a) ji + n~ pairs, b) e + e~ 
pairs. Both particles are detected in the central region. The full lines are the results of a fit using 
a Gaussian distribution for the signal region (convoluted with an exponential tail to account 
for energy loss in the case of di-electron decays) and an exponential distribution for the non- 
resonant background. The mass spectra from two photon processes (LPAIR simulation) are 
shown as hatched histograms. is the number of J ftp events obtained from the fit. 
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Figure 3: Mass spectra for events of the J /if) selection in the low W analysis: a) pairs 
for the FMD-FMD sample and b) for the FMD-CTD sample. The open histogram repre- 
sents the prediction of a Monte Carlo simulation including diffractively produced J /if) mesons 
(DIFFVM) and muon pairs from two photon processes (LPAIR) which are also shown sepa- 
rately as hatched histograms. 




Figure 4: Control distributions for the quasi-elastic J /if) selection (data set I. in Table p. a) 
Polar angle and b) energy of the scattered positron, c) polar angle and d) momentum of the 
reconstructed J /if) candidates. The error bars on the data points are statistical only. Results of a 
diffractive Monte Carlo simulation (DIFFVM) normalised to the data are shown as histograms. 
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background under the J /if) peak is determined by fitting the sidebands using an exponential 
distribution and is found to be 12 ± 3% on average. The error includes the uncertainties of the 
resonance parameterization and of the background shape where the latter was estimated using a 
power law as alternative. The non-resonant background is mainly due to dilepton production by 
two photons as simulated by the Monte Carlo generator LPAIR, but there is also a contribution 
from hadrons misidentified as leptons. Some distributions for events in the signal region are 
shown in Fig. Q as well as the predictions from a diffractive simulation, DIFFVM. The dis- 
tributions for the scattered positron and for the reconstructed J ftp meson are reasonably well 
described by the simulation. Remaining differences between data and simulation were checked 
to have a small impact on the results and are accounted for by the systematic uncertainty. 

The distribution of the quasi-elastic J/ip events from the three data sets (CTD-CTD, FMD-CTD 
and FMD-FMD) in the kinematic plane x versus Q 2 is displayed in Fig. |5] before applying Q 2 
and W cuts. 

4.1 Cross Sections as Functions of W and Q 2 

In order to measure the elastic cross section, the quasi-elastic data sample is divided into two non 
overlapping classes, forward untagged and forward tagged (see Table [IJ, I.), which are enriched 



m J/f -> ee (CTD-CTD) A J/V v fJ,/J. (FMD-CTD) 
O J/ty > fifi (CTD-CTD) ★ J/i> fj,[L (FMD-FMD) 




1 10 10 

Q 2 [GeV 2 ] 

Figure 5: Distribution of quasi-elastic J/ip candidates in the kinematic (x, Q 2 ) plane (data sets 
I. and II. in Table [1]), before applying Q 2 and W cuts. Lines of constant W, energy E e and polar 
angle 6 e of the scattered positron are shown and the different analysis regions are indicated. 
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in elastic and proton dissociation processes, respectively. The data are binned in Q 2 and W and 
the non-resonant background is determined and subtracted for each bin. The "true" numbers 
of elastic and proton dissociative events N ei and N pd are extracted by unfolding them from the 
number of events with and without a tag of the forward detectors (see section [I]). The efficiencies 
for tagging and for non-tagging of elastic and proton dissociative events are determined by 
studying the response of the forward detectors and are incorporated in the detector simulation. 
The tagging efficiency for proton dissociative events with M x £1.6 GeV is found to be 92% 
on average. Note that in this procedure no assumption is made for the absolute or relative cross 
sections of the two processes. 

A small correction (~ 3%) due to the presence of spurious hits in the FMD which are not 
described by the Monte Carlo simulation is applied. Further corrections account for the con- 
tamination from decays of the ip(2S) meson into a J/ij; and undetected neutral particles (based 
on the results of section |5|) and for initial state radiation (section |37J). Using the integrated 
luminosity and the sum of the branching fractions for the J/ip meson to decay into ji + ^ or 
e + e~, an integrated ep cross section is calculated for each Q 2 and W bin. 

In the Born approximation, the electroproduction cross section is related to the 7*p cross section 
by 

dMe \y^ v) = Mt*p - Wp) = rv T ( 7 *p - jam(i + sR), (7) 

where R = a L /a T , o T and o L are the transverse and longitudinal 7*p cross sections. Y is the 



flux of transverse virtual photons [ p9| ] and e is the flux ratio of longitudinally to transversely 
polarized photons, given by 

r = ^7^-(i + (i-y) 2 ); e= (8) 

2iryQ 2 l-y + y 2 /2 

Virtual photon-proton cross sections are computed using equation (0) after integrating over the 
Q 2 and W bins used in the analysis. The difference between a(YP ~^ J/i>p) — °t + and 
Vtot{l*P — > J/4>p) = cr T + a L is negligible here since (s) = 0.99. 

For the analysis in which both tracks are detected in the central detector the systematic uncer- 
tainties of the cross sections are estimated to be 17% in total, and are only slightly dependent 
on the kinematics. They consist of uncertainties due to detector efficiencies and resolution 
(10%), uncertainties in the estimation of background (11%, dominated by proton dissociation 
and ip(2S) decays), radiative corrections and bin centre determination (4%), the J/ip decay 
branching ratio, and luminosity determination (4%). Part of the systematic error (9%) affects 
only the overall normalization. The uncertainty arising from the proton dissociation background 
is estimated by varying the cuts to select proton dissociation events, by changing the M x de- 
pendence assumed in the Monte Carlo simulation, and by changing the model used for the 
fragmentation of the system X. The uncertainty due to the subtraction of non-resonant back- 
ground is determined by varying the assumed shape of the background and using alternative 
methods for its determination, such as sideband subtraction. 

For the low W analysis (one muon in the FMD) a different procedure to extract the cross section 
was adopted due to limited statistics. Since the contribution of hadrons misidentified as muons 
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50 100 

W [GeV] 

Figure 6: Cross sections for elastic J / ip production as a function of W at different values of Q 2 , 
measured at HERA in this and other analyses [0, ^ |], |5|] . Data for Q 2 > have been scaled by 
factors 5, 50 and 100 as indicated. The inner error bars on the points from this analysis indicate 
the statistical errors, while the outer bars show the statistical and systematic uncertainties added 
in quadrature. The dash-dotted lines are the results of fits of the form W s to the data for each 
Q 2 . For Q 2 > the fits are for HI data only. The full and dashed lines are results of calculations 
from Frankfurt et al. [§] using different parameterizations of the proton gluon densities. The 
small arrows at W = 200 GeV indicate the sensitivity of this prediction to a change of the 
charm quark mass from 1.4 GeV to 1.5 GeV. 

is negligible here, the non-resonant background is subtracted using the LPAIR Monte Carlo 
simulation. Alternatively it is estimated from the sidebands of the mass spectrum. The proton 
dissociation background is subtracted assuming the same fraction as determined in the CTD- 
CTD analysis; this assumption was verified by comparing the response of the forward detectors. 
The efficiency of the FMD is determined using a sample of J/ij) photoproduction events and 
is cross checked by the control sample with both muons in the FMD. On average, the FMD 
efficiency is found to be ~ 81%. The total systematic uncertainty of the cross section in the 
low W analysis is 25%, dominated by the uncertainties in the subtraction of non-resonant and 
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proton dissociation backgrounds and by the uncertainty of the FMD efficiency. 

The 7*p cross sections are shown in Fig. ^| and are given in Table |] as functions of W in three 
bins of Q 2 (2 < Q 2 < 6 GeV 2 , 6 < Q 2 < 18 GeV 2 and 18 < Q 2 < 80 GeV 2 ). Also shown 
are measurements of the ZEUS collaboration^ at similar values of Q 2 . The cross sections of the 
present analysis are quoted at values of W and Q 2 after applying bin centre corrections using 
the measured W and Q 2 dependences. The W dependence, which in pQCD based models is 
related to the x dependence of the gluon density in the proton, is found to be similar to that 
obtained in the photoproduction limit at HERA (also shown in Fig. When parameterized in 
the form W 5 , fits to HI and ZEUS photoproduction data yield 8 = 0.77 ± 0.18, while the HI 
data for Q 2 > yield 5 = 0.84 ±0.20 at Q 2 = 3.5 GeV 2 and 5 = 1.3 ±0.4 at Q 2 = 10.1 GeV 2 , 
where the errors include statistical and systematic uncertainties. The fits are shown in Fig. |£]. 

In Fig. ^ also predictions of the model by Frankfurt et al. |Q] are included. Gluon densities 
from GRV(HO) |Q and MRSR2 [@4]] are used at an effective scale depending on Q 2 and on 
the separation of the quarks within the J/ip. The prediction using MRSR2 describes the slope 
of the data well while the calculation using GRV(HO) is too steep at low values of Q 2 . At 
small Q 2 , the absolute magnitudes of the predictions are very sensitive to the input value for 
the charm quark mass (m c = 1.4 GeV was chosen here), as is indicated by the arrows in Fig. |6|: 
For Q 2 = and W = 200 GeV, for example, a change from m c = 1.4 GeV to m c = 1.5 GeV 
reduces the prediction by more than 40%. 

The Q 2 dependence of the cross section for W = 90 GeV (40 < W < 160 GeV) is shown in 
Fig. [7] and given in Table |5|. It is well described by a fit oc (Q 2 ± m^) - " with n = 2.38 ± 0.11. 
In order to study a possible change in the observed Q 2 dependence, which may be an indication 
for the importance of non-perturbative effects, the fits are repeated in two Q 2 regions leading to 
n = 2.12 ± 0.20 for Q 2 < 12 GeV 2 and n = 2.97 ± 0.51 for Q 2 > 12 GeV 2 . The errors contain 
statistical and systematic uncertainties. In Fig. [7] the model of Frankfurt et al. with the MRSR2 
gluon distribution which was seen to give a good description of the W dependence (Fig. |5|) is 
also compared to the data. The Q 2 dependence is reasonably well described by the prediction. 

4.2 t Distribution and Elastic Slope Parameter 

The elastic slope parameter b is determined assuming that the t dependence of the elastic J/ip 
cross section can be parameterized by a single exponential distribution e bt . Three contributions 
are fitted to the forward untagged J/tp sample, corrected for acceptance, losses and smearing 
effects. These are: 

• One exponential distribution e bt with a slope b as free parameter describing elastic J / ip 
production. 

• The non-resonant background is described by the sum of two exponential distributions 
contributing in total 12%. The t- slopes of the non-resonant background are determined 

5 In the present paper, we do not use fixed target data for comparison because experimental conditions and 
methods are different and lead to uncertainties in the comparison: for example most experiments used heavy 
nuclei as targets and define elastic processes differently. 
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Figure 7: The cross section for elastic J/ift production at W = 90 GeV as a function of Q 2 . The 
inner error bars on the points from this analysis indicate the statistical errors, while the outer bars 
show the statistical and systematic uncertainties added in quadrature. Also shown are previous 
measurements in photoproduction (indicated by Q 2 = 0) [Q, |5p and deep inelastic scattering 
[0 0]. The full line is a fit of the form (Q 2 + mV)~ n , yielding the result n = 2.38 ± 0.11. The 
dashed line is the prediction of Frankfurt et al. [jjjl using the MRSR2 [^TJ] gluon density. 

using the sidebands of the J/ip mass distribution. The background fraction depends 
strongly on t. For the estimation of the systematic uncertainty, the total amount is varied 
within the range 5 — 16%. 

• The proton dissociation background is described by one exponential with a slope param- 
eter 1.4 GeV -2 . This is compatible with studies of the forward tagged data set taking into 
account non-resonant background. The total contribution is fixed to 13% and is varied be- 
tween 5% and 28% to estimate the systematic error, while the slope was varied between 
0.8 GeV -2 and 2.0 GeV -2 . 

No correction for background from ip(2S) decays is applied since the total contribution is small 
at low Q 2 . The result of the fit which is carried out up to \t\ = 1.2 GeV 2 is shown in Fig. || 
(X 2 /NDF = 1.6/4). The elastic slope parameter is 

b = 4.1 ± 0.3 (stat.) ± 0.4 (syst.) GeV~ 2 (9) 
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Figure 8: \t\ distribution for the forward untagged J/ip sample, corrected for acceptance, losses 
and smearing effects. The dashed line is the result of a fit taking the background contributions 
into account as described in the text. The full line corresponds to the elastic contribution assum- 
ing an exponential distribution. The contributions from proton dissociation and non-resonant 
background are shown separately. The error bars on the data points are statistical only. 

for mean values (W) = 96 GeV and (Q 2 ) = 8 GeV 2 . The systematic uncertainty was estimated 
by varying the fit range by ±0.4 GeV 2 and by varying the background contributions and the 
corresponding slopes within the ranges given above. 

This result for b is compatible with the values obtained by HI [Q| and ZEUS [£]] for elastic J/ip 
photoproduction at similar values of W: b = 4.4 ± 0.3 GeV 2 (HI) and b = 4.6 ± 0.6 GeV 2 
(ZEUS), as well as the ZEUS measurement [[| of b = 5.1 ± 1.3 GeV 2 for 2 < Q 2 < 40 GeV 2 
and 55 < W < 125 GeV. With the present statistics no significant dependence of the b- 
parameter on W is found (see Table |2|); there is however an indication for a decrease of b with 

Q 2 . 

4.3 Decay Angular Distributions for Quasi-Elastic J/ 1\) Production 

In order to investigate the helicity structure of J/ip meson production the angular distribu- 
tions of the decay leptons in the helicity frame are used. In this frame, the J/ip direction in the 
7*p centre of mass system serves as the quantisation axis. Three angles are defined: the polar 
(9*) and azimuthal {ip) angles of the positive decay lepton in the J/ip rest frame. The third angle 
is the angle <p between the normals to the J/ip production plane (defined by the J/ip and the 
scattered proton) and the electron scattering plane in the 7*p centre of mass system. 
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The one-dimensional distributions in cos 9* and the polarization angle \1> = ip — <p are extracted. 
If the helicity of the virtual photon is retained by the J/ip meson (s-channel helicity conserva- 
tion hypothesis, SCHC), the full angular distribution is a function of cos 9* and \P only. 

The acceptance corrected cos 9* and \& distributions are shown in Fig. Both the forward 
tagged and untagged event samples are used, since the helicity structure is expected to be inde- 
pendent of whether the proton dissociates or not, a hypothesis supported by the data. The cos 9* 



40 < W < 
2 < Q 2 < 8 GeV 2 


160 GeV 

8 < Q 2 < 80 GeV 2 


b = 4.4 ± 0.4 ± 0.4 GeV" 2 


b = 2.5 ± 0.6 ± 0.6 GeV" 2 


2 < Q 2 < 
40 < W < 100 GeV 


80 GeV 2 

100 < W < 160 GeV 


b = 4.0 ± 0.4 ± 0.4 GeV" 2 


b = 4.1 ±0.5 ±0.5 GeV" 2 



Table 2: Slope parameters b of the elastic J/ip meson t distribution for different Q 2 and W 
domains. 




Figure 9: Angular distributions for the positive J/tp decay lepton in quasi-elastic production 
processes e + p -> e + J/ip + X at 40 < W < 160 GeV. a) cos 9* for 2 < Q 2 < 8 GeV 2 and 
b) for 8 < Q 2 < 80 GeV 2 ; c) and d) the polarization angle ^ in the same Q 2 regions. The inner 
error bars indicate the statistical errors, while the outer bars show the statistical and systematic 
uncertainties added in quadrature. The lines are fits to the data as described in the text. 
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distribution is related to the spin density matrix element r^, the probability of the J/ip meson 
to be longitudinally polarized, according to 

<1<T ocl+r 04 + (l-3r° 4 )cos 2 r. (10) 



d cos 9* 

A fit to the data shown in Figs. ^| a) and b) yields 

r° 4 = 0.15 ±0.11 for (Q 2 ) = 4 GeV 2 , (11) 

r 04 

'oo 



04 - 0.48 ±0.15 for {Q 2 ) = 16 GeV 2 . (12) 



The ^ distribution is related to the spin density matrix element rj_ x : 

da 



OC 



1 -er 1 1 _ 1 cos2^. (13) 



A fit to the data shown in Figs. ^|c) and d) yields 

r 1 1 _ 1 = 0.50 ±0.08 for (Q 2 ) = 4 GeV 2 , (14) 
r 1 1 _ 1 = 0.39 ±0.13 for (Q 2 ) = 16 GeV 2 . (15) 

In the case of s-channel helicity conservation and natural parity exchange (NPE) the matrix 
elements r° 4 and r{_ x are related by r 4 _ x = |(1 — r° 4 ). Using this relation and the measured 
values for Tq 4 , one obtains values for r\_ x which agree to within one standard deviation with 
those obtained from the ^ angular distributions, thus supporting the SCHC and NPE hypothe- 
ses. 

Under the assumption of SCHC the measurement of the r[} 4 matrix element can be used for the 
determination of R, the ratio of the longitudinal to the transverse cross section: 



a T el- rgo 



R=~ = L T^M- (16) 



Using this relationship and the cos 9* distribution R is determined in two Q regions: 

R = 0.181^ for (Q 2 )=4GeV 2 , (17) 
R = OMtoll for (g 2 ) = 16GeV 2 . (18) 

A measurement of R = 0.41+^;^ by the ZEUS experiment [§] at (Q 2 ) = 5.9 GeV 2 and 
(W) = 97 GeV is compatible with these values. Taking into account the photoproduction 
measurements of R = 0.17 ± 0.14 [Q] and R = -0.01 ± 0.09 [§], which are compatible with 
the expectation R = for Q 2 = 0, a rise of R with increasing Q 2 is suggested by the data. 

The measured values of the R parameter are significantly smaller for J/ip than for elastic p 
meson production at HERA [Q, ^ Q at similar Q 2 ; but they are of the same order if compared 
at the same value of Q 2 /m 2 /, where my is the mass of the p or the J/ip. 



22 



5 Quasi-elastic if){2S) Production 



For the selection of ip(2S) mesons the decay channel ip(2S) — > J/ipn + n~, where the J /if) 
decays either in two electrons or two muons, is used. In this case no separation between elastic 
and proton dissociation is attempted due to limited statistics. The goal is to derive the ratio of 
cross sections for J /if) and ip(2S) production as a function of Q 2 . The lower Q 2 cut is reduced 
to 1 GeV 2 since the Q 2 dependent acceptance corrections cancel almost completely in the cross 
section ratio. 

The signals in the quasi-elastic if)(2S) selection are displayed in Fig. [T0|. For the determination 
of the ip(2S) to J /if) ratio the non-resonant background is subtracted using the sidebands of the 
di-lepton mass spectrum in the case of the J/tp meson and using the sidebands of the Am = 
m 4>(2S) ~ m 4> distribution for the if>(2S) meson. The data are divided in three Q 2 bins: 1 < 
Q 2 < 5 GeV 2 , 5 < Q 2 < 12 GeV 2 and 12 < Q 2 < 80 GeV 2 . The cross section ratio is shown 
in Fig. [TT]. The systematic uncertainty of the ratio amounts to 16% in total and is dominated by 
the contribution from the track reconstruction efficiency for the low momentum 7r + 7r pair. 



The measurement at low Q 2 agrees well with the HI photoproduction measurement p4|]. An 
indication of a rise with Q 2 at the level of two standard deviations is observed which is also 
predicted in models by Frankfurt et al. [^] and Nemchik et al. [|TT)]. In ref. [§] an asymptotic 
value of (t^(2S)/<J^ « 0.5 is expected for Q 2 3> m|. 



6 Inclusive and Inelastic J/ij) Production 

Inclusive J /if; production is studied in the kinematic range 2 < Q 2 < 80 GeV 2 and 40 < W < 
180 GeV covering 0.2 < z $ 1.0 for the muonic decay of the J /if), while for the decay to 
electrons z > 0.5 is required. The restricted z region for J /if) — > e + e~ is due to the smaller 
acceptance for electrons and larger background at low z values. The elasticity z (defined in 
section [O) denotes the ratio of energies of the J / if> and of the exchanged photon in the proton 



rest frame. Two sets of differential cross sections are determined. First an inclusive cross 
section is derived where in the given kinematic region all J / if> mesons are selected irrespective 
of the production mechanism, thus including inelastic and elastic contributions. The inclusive 
cross sections are compared to the predictions of the Soft Colour Interaction Model [EUR. A 



second set of differential cross sections is derived for inelastic J / if) production which can be 
compared to the predictions within the NRQCD factorization approach [ |19| ] containing colour 
octet contributions. 

An "inelastic" production process can be defined experimentally in several ways. In previous 



photoproduction analyses [Q, [14]] cuts in the variable z, e.g. z £ 0.9, were used to suppress 
elastic and proton dissociative events. In the present analysis a different approach is chosen 
because colour octet contributions are, in leading order in a s , predominantly expected at large 
z. This is because the cc pair can be produced with no other particles in the final state, i.e. z ~ 1 
(see Fig. |I]e). Its non-perturbative evolution into the J /if) meson reduces the value of z only 
slightly, and applying a z cut as done previously would reduce the expected colour octet con- 
tributions together with the quasi-elastic ones by an unknown amount. In the present analysis a 
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Figure 10: a) Mass spectrum of the four particles £ + 1~tt + it~ and b) Am = m^s) 
the ip{2S) candidate events, i.e. for events with \mu — m$\ < 300 MeV. The DIFFVM Monte 
Carlo simulations for the signals are shown for comparison. 
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Figure 1 1 : Ratio of cross sections for the quasi-elastic production of ij)(2S) and J / tp mesons as 
a function of Q 2 for this analysis and for the HI photoproduction measurement [JPJ], corrected 
for the most recent branching fraction BR(ip(2S) -> J/ipn+n-) = 30.2 ± 1.9% @. The 
inner error bars on the HI points from this analysis indicate the statistical uncertainty, while the 
outer bars show the statistical and systematic uncertainties added in quadrature. The prediction 



from based on colour dipole phenomenology is also displayed. 
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a) b) 
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Figure 12: Di-lepton mass spectra for events of the inclusive (points) and inelastic (histogram) 
J/ip selection: a) 0.2 < z < 0.6, b) 0.6 < z £ 1. The curves are the results of fits of Gaussian 
distributions for the signal (convoluted with an exponential tail to account for energy loss in the 
case of di-electron decays) and an exponential distribution for the non-resonant background. 



cross section is determined suppressing contributions of low mass Mx following a suggestion 
of JT9|] . This suppression of low masses is achieved indirectly by requiring a minimal calori- 
metric energy in the forward region of the detector. This requirement selects high masses and 
suppresses elastic and proton dissociative events characterised by small M x corresponding to 
small energy deposits in the forward direction. Colour octet contributions are expected to have 
M x % 15 GeV M and are retained. 



6.1 Data Analysis 

The selection criteria as described in Table |T] are used. The di-lepton mass spectra of the selected 
events are shown in Fig. [12| separately for 0.2 < z < 0.6 and 0.6 < z ^ 1, both for the inclusive 
and the inelastic selection (Ef w d > 5 GeV). Since the non-resonant background increases with 
decreasing z the background fraction is determined from the mass spectra in bins of z by fitting 
the signal and background as in section^. For the determination of the differential cross sections 
a correction is applied according to the z values of the events. 

Acceptance and efficiencies are determined using a simulation tailored to describe the data, 
which consists of a mix of diffractive events generated by DIFFVM [R3p and of inelastic events 



generated by EPJPSI [p6[| . The diffractive events are composed of elastic and proton disso- 
ciative contributions in a ratio consistent with the signals observed in the forward detectors 
(compare section |j). EPJPSI generates events according to the Colour Singlet Model. Both 
models were previously shown to describe quasi-elastic and inelastic data respectively (see for 
example Contributions from other processes such as 6-decays or the hadronic component 



of the photon are expected to contribute only for z ^ 0.4 and are estimated to be negligible. 

The EPJPSI contribution is normalized to the data in the interval 0.4 < z < 0.8 and the 
DIFFVM contribution is added to describe the data in the region z > 0.95 (compare also 
Fig. p3|d). Numerous checks were carried out to ensure that all important aspects of the data are 
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Figure 13: Comparison between data and Monte Carlo simulations for inclusive J / ip production 
after all selection cuts in Table |T] and after background subtraction. Shown are distributions of 
a) Q 2 , b) W, c) the square of the J/ij) transverse momentum in the laboratory frame pf^,d) the 
elasticity z, e) the energy Ej w d deposited in the LAr calorimeter for 9 < 20°, and f) the energy 
E CO ne in a cone with radius R = a/ (Ar/) 2 + (A0) 2 = 1 (rj = — lntan(6 l /2)) around the J ftp 
direction of flight. The results of the combined Monte Carlo simulation (DIFFVM and EPJPSI, 
full lines) and of the DIFFVM simulation only (dashed lines) are shown. The error bars on the 
data points are statistical only. 
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Figure 14: Distribution of simulated J/ip events as a function of the generated value of the mass 
M x . The mixed Monte Carlo sample (DIFFVM and EPJPSI) is shown before (full histogram) 
and after (dashed histogram) applying the cut Ef wd > 5 GeV. The diffractive contribution as 
simulated by DIFFVM before the cut is also shown (dotted histogram). 



well described by this mix. Comparisons between the data and the Monte Carlo simulation are 
shown in Fig. |I3| . 

The systematic errors in this analysis are typically 18% and are dominated by uncertainties in the 
acceptance corrections, mainly due to the model uncertainty of the Monte Carlo description, the 
subtraction of non-resonant background, and reconstruction efficiencies. The largest systematic 
uncertainty (up to 32%) is found for large z and small p 2 ^ values. 

In a second step inelastic cross sections are determined for events with a large energy deposition 
in the forward region of the LAr calorimeter, namely Ej wd > 5 GeV for polar angles 9 < 20°. 
This requirement is an indirect cut on the mass of the hadronic system X. Its effect can be seen 



for simulated events in Fig. [14] where the M x distribution is shown for the mixed simulation 
(EPJPSI+DIFFVM) and for the fraction of the diffractive simulation (DIFFVM) separately. 
The latter dominates at low values of M x and is suppressed efficiently by the cut on Ef w d. 
The differential cross sections for M x > 10 GeV are thus determined by applying the cut 
Ef w d > 5 GeV and then correcting to M x > 10 GeV using the Monte Carlo simulation. 



6.2 Differential Cross Sections 

Inclusive Cross Sections and Soft Colour Interactions Differential ep cross sections for 
inclusive J / ip production are given in Table ^ and shown in Fig. [15] (open points) as functions of 
Q 2 , p 2 ^, z, y* (the rapidity of the J/ip in the j*p centre of mass system) and W. The prediction 
of the Soft Colour Interaction Model (dotted histogram in Fig. [15]) which is computed using a 
modified version of AROMA [ETL f|5J][], is compared to the data. Although the model gives a 



6 The following parameters are used in addition to standard settings: charm mass m c = 1.4 GeV, GRV(HO) 
parton densities and R = 0.5, where R parameterizes the probability of a colour exchange between partons. 
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Figure 15: Differential cross sections for the inclusive (open points) and inelastic (M x > 
10 GeV, black points) ep — > e J/ip X process, a) da/dQ 2 , b) da/dp 2 ^ (see also footnote 
H concerning the theoretical prediction), c) da/dz, d) da/dy* and e) da/dW. The kinematic 
region is 2 < Q 2 < 80 GeV 2 , 40 < < 180 GeV and 2 > 0.2. The inner error bars are 
statistical, the outer error bars contain statistical and systematic uncertainties added in quadra- 
ture. The dotted histogram gives the prediction from the SCI model in AROMA 2.2 [ ]2~T| , p3| ] for 
inclusive J/ip production. The curves are predictions for inelastic J ftp production within the 



NRQCD factorization approach [|19| for the colour singlet contribution (dashed line) and the 
sum of singlet and octet contributions (full line). 
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reasonable description of the shapes of several distributions, there are major discrepancies in 
the z distribution and in the absolute values of the measured and predicted cross sections. For 
small z, the SCI model as implemented in AROMA is expected to fall below the data due to 
the missing hard contributions of the Colour Singlet Model which should also be taken into 
account [|3p. At large z the AROMA SCI prediction is below the measured cross section by 
approximately a factor of four. 



Inelastic Cross Sections and NRQCD Factorization Approach The differential cross sec- 
tions for inelastic J/ij) production, that is for Mx > 10 GeV, are also displayed in Fig. [L5] (full 
points). In comparison with the inclusive cross sections the effect of requiring a high mass is 
most significant in the shapes of the distributions of z and y*. 



The results of the calculations by Fleming and Mehen [|T9p who applied the NRQCD factoriza- 
tion approach to electroproduction of J/ij) mesons are shown in Fig. [15] for comparison. The 
predicted cross sections include the contributions from the colour octet states 3 P , 1 S which 
are of order 0(a s ) and the colour singlet state 3 S'i (of order 0(a^))[]. The sum of these contri- 
butions shown in Fig. [T5J is computed using GRV(LO) [ fTO| ] parton densities; the colour singlet 
contribution is also shown separately. Note that these predictions are for cc pairs from the hard 
subprocess and do not include any hadronisation effects. The hadronisation of the colour octet 
cc pairs into a colour singlet J ftp is believed to proceed via emission of soft gluons 0. 

The colour octet contribution dominates the cross section for all Q 2 (Fig. p3|a). The colour 



singlet contribution (dashed curves in Fig. |T5|) is seen to fall below the data by factors 2—3 while 
the prediction for the sum (full curves in Fig. [T5|) is overall too large in absolute magnitude by 
up to a factor 3. The shapes of the data distributions are not well reproduced by the calculation: 
the predicted Q 2 and p 2 ^ distributions are steeper than the data and the y* distribution increases 
towards larger values of y* instead of falling. The W distribution agrees in shape but overshoots 
the data. The full prediction of the z distribution is at present not calculable [[T^] and is therefore 
not shown. 

The observed differences in magnitude between the predicted and measured cross sections prob- 
ably call for an overall adjustment of the fitted transition matrix elements while the shapes may 
be influenced by a relative adjustment of the individual contributions. There is a hint that these 



differences increase towards low Q (see Fig. |T5|a, full points and full curve). The theoretical 



predictions are also expected to be more precise for larger Q 2 and p 2 ^ [|TP|]. Therefore the com- 
parison was repeated for Q 2 > 4 GeV 2 and p 2 ^ > 4 GeV 2 (see Table [7]), but no significant 
change in the conclusions was found. 

7 Spectroscopic notation is used: 2S+1 Lj where S, L and J denote spin, orbital and total angular momentum of 
the cc system that is produced in the hard process. The following values for non-perturbative long range transition 
matrix elements are used: (O^ ( 3 Sx)) = 1.1 GeV 3 , (O^^Sq)) = 0.01 GeV 3 , and (O J { ^( 3 P ))/m 2 c = 

0.005 GeV 3 . The octet matrix elements (O^) were estimated from fits to the CDF data performed in [^6| while 
the singlet matrix element is calculated from the measured electronic decay width of the J/ip. 

',■>!> 



8 The sharp edge observed in Fig. |l5| b in the theoretical p 2 , curve is a consequence of the missing hadronisa- 



tion. 
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Data set 


a(ep — > e J/% 
Q 2 > 2.0 GeV 2 


b X) [nb] 

Q 2 ,^>4.0GeV 2 


Inclusive 

M x > 10 GeV 
z < 0.9 


1.30 ±0.06 ±0.24 
0.51 ±0.04 ±0.09 
0.46 ±0.04 ±0.09 


0.50 ±0.04 ±0.09 
0.17 ±0.02 ±0.03 
0.13 ±0.02 ±0.02 



Table 3: Integrated cross sections for 40 < W < 180 GeV, Q 2 < 80 GeV 2 and z > 0.2 in 
two kinematic regions, Q 2 > 2.0 GeV 2 and Q 2 and p 2 ^ > 4.0 GeV 2 . Results are given for the 
inclusive selection and for two inelastic selections, M x > 10 GeV and z < 0.9. 



6.3 Integrated Cross Sections and Comparison with Photoproduction 

In Table [5] the integrated cross sections for e + p — > e + J ftp + X in the kinematic region 
2 < Q 2 < 80 GeV 2 , 40 < W < 180 GeV and z > 0.2 are summarised. They are given for the 
inclusive selection, for the inelastic selection corresponding to Mx > 10 GeV, and for z < 0.9 
as in previous photoproduction analyses [Q, [T3j]. In addition, the cross sections after imposing 
the additional cuts Q 2 > 4.0 GeV 2 and p\^ > 4.0 GeV 2 are provided. 

The total cross section for 7*p — ► J/ip X is computed according to equation (0) as a function 
of the 7*p centre of mass energy W at (Q 2 ) = 9 GeV 2 and is given in Fig. [16] and Table [8|. The 
cross section is determined for the inclusive data 0.2 < z £ 1 and, in view of a comparison 
with photoproduction, also for an inelastic selection using a cut z < 0.9 f\. The W dependence 



is seen to be very similar to that in the photoproduction data [Q, [140 also shown in Fig. |16 
The W dependence, parameterized as W 5 , yields 5 = 0.95 ±0.11 for the inclusive data and 
5 = 0.89±0.20 for the data with z < 0.9. The photoproduction data, including HI and ZEUS, 
are described by 5 = 0.91 ± 0.26. The errors include statistical and systematic uncertainties. 



6.4 Decay Angular Distributions 

Measuring the polarization of the J/ip is thought to be a way of distinguishing the various 
contributions to J / ip production. The polar (9*) decay angular distributions in the helicity frame 
are shown in Fig. [T7]for the fully inclusive case and for the inelastic selection M x > 10 GeV, 
in the kinematic region 2 < Q 2 < 80 GeV 2 , 40 < W < 180 GeV and z > 0.2. 

The cos 9* distribution is predicted to have the form 

, da - oc 1 + a cos 2 fl*. (19) 
a costr 

For J/ij) production via the colour singlet mechanism a ~ 0.5 is expected for the kinematic 



range studied here [[47]]. If colour octet contributions are present, |a| £ 0.5 is expected, where 
a can be negative, zero or positive depending on which intermediate cc state dominates the 
production [|I9]]. 

The data yield values of a = 0.541^26 for the inclusive case and a = 0.771° 38 f° r tne inelastic 
selection (Mx > 10 GeV), including statistical and systematic uncertainties. The uncertainties 
are too large to draw definite conclusions. 

9 The photoproduction data are given for < z < 0.9. This was achieved by an extrapolation to z = assuming 
contributions from photon gluon fusion only. This contribution at small z is however negligible in the comparison. 
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Figure 16: Total cross sections for 7*p — > J/tp X from this analysis at (Q 2 ) = 9 GeV 2 . The 
inclusive cross section (0.2 < z £ 1.0) is shown as a function of W (multiplied by a factor 
1.5 for clarity), as well as the cross section for 0.2 < z < 0.9. Photoproduction data [|J [TJ] 
with similar cuts in z are included for comparison. The inner error bars on the points from 
this analysis indicate the statistical uncertainty, while the outer bars show the statistical and 
systematic uncertainties added in quadrature. 
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Figure 17: Differential cross sections da /d cos 9* for ep — > e X in the kinematic region 
2 < Q 2 < 80 GeV 2 , 40 < < 180 GeV and z > 0.2. The inclusive cross section and the 
inelastic cross section (Mx > 10 GeV) are shown. The inner error bars indicate the statistical 
uncertainty, while the outer bars include the statistical and systematic uncertainties added in 
quadrature. The lines are fits to the form ~ 1 + a cos 2 9*. 
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7 Summary and Conclusions 



Measurements of elastic J/tp production in deep inelastic scattering (DIS) with 25 < W < 
160 GeV and 2 < Q 2 < 80 GeV 2 have been presented. They are more precise and cover a 
larger kinematic range than previous analyses at HERA. The dependence of the cross section 
a(Y? J /i> p) onW is found to be proportional to W 5 with S ~ 1, as was also observed in 
photoproduction. The Q 2 dependence is measured to be oc l/(Q 2 + m 2 b ) n withn = 2.38 ±0.11. 
Both the W and Q 2 dependence are well described by a model based on perturbative QCD. 

Assuming that the t dependence of elastic J/ip production can be described by one exponential 
distribution, the slope parameter is determined to be b = 4.1 ± 0.3 (stat.) ± 0.4 (syst.) GeV -2 , 
compatible with the value found in photoproduction. The helicity structure of quasi-elastic 
J/ijj production in DIS has been investigated and no evidence is found for a violation of s- 
channel helicity conservation. Assuming SCHC the ratio R of the longitudinal to the transverse 
cross section has been determined using the cos#* distribution in two Q 2 regimes; the result is 
R = 0.181^ for (Q 2 ) = 4 GeV 2 and R = 0.94l°;^ for (Q 2 ) = 16 GeV 2 , suggesting a rise 
with Q 2 . 

The first evidence from HERA for quasi-elastic production of ip{2S) mesons in DIS has been 
reported. The increase of the ratio of cross sections for ip(2S) and J/ip production with Q 2 
predicted by models is supported by the data. 

Data have been presented for the inclusive production of J/ip mesons in deep inelastic scatter- 
ing, covering the kinematic region 40 < W < 180 GeV, 2 < Q 2 < 80 GeV 2 and 0.2 < z < 1. 
Differential ep cross sections are computed as functions of Q 2 , p 2 ^, z, y* and W. The model 
of Soft Colour Interactions, a non-perturbative phenomenological approach to the description 
of inclusive J/ip production, is compared to the data. The dependences of the differential cross 
sections on several variables are reasonably well described by the model, but the normalisations 
and the z dependence are not reproduced. 

Using a selection cut designed to reject events with a low mass hadronic system, diffractive 
events are suppressed and inelastic cross sections (M x > 10 GeV) are extracted. A lead- 
ing order calculation in the NRQCD factorization approach using long range matrix elements 
determined from J / ip production in pp collisions at the Tevatron is confronted with our mea- 
surements of differential inelastic cross sections. The shape and magnitude of the differential 
distributions are not described by the theoretical prediction. The comparisons may indicate the 
need to decrease the size of the colour octet long distance matrix elements or to change the 
relative importance of the different colour octet contributions, and/or to include higher orders 
in the NRQCD perturbative expansion. The colour singlet contribution alone is below the data 
by factors 2 — 3. 
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in 2 \ 


W interval 


\ w ) 




a(j*p 




[GeV 2 l 
IV C v J 




V c J 


fGeVl 
IV c v J 






[nb] 




25 


- 40 


32.0 


11.7 


± 3.2 (stat.) ± 2.9 (syst.) 




40 


- 60 


49.3 


22.5 


±4.0 


±3.6 


3.5 


60 


- 80 


69.5 


26.3 


±3.9 


±4.2 


80 


- 100 


89.6 


33.1 


±5.1 


±5.2 




100 


- 120 


109.6 


30.7 


±5.8 


±4.8 




120 


- 160 


138.6 


54.9 


±8.9 


±8.8 




40 


- 80 


57.5 


5.4 


± 1.3 


±0.9 


10.1 


80 


- 120 


98.4 


10.3 


±2.2 


± 1.7 




120 


- 160 


138.6 


17.8 


±4.2 


±3.1 


33.6 


40 


160 


84.4 


1.34 


±0.37 


±0.24 



Table 4: Cross sections for the elastic process 7*p — > J/ip p in bins of W for three Q 2 regions: 
2 < Q 2 < 6 GeV 2 , 6 < Q 2 < 18 GeV 2 and 18 < Q 2 < 80 GeV 2 . 



Q 2 interval 


(Q 2 ) 


cr(7*p 


-> J/rpp) (W 


= 90 GeV) 


[GeV 2 ] 


[GeV 2 ] 




[nb] 




2.0 - 3.2 


2.6 


31.9 


± 2.5 (stat.) 


± 5.1 (syst.) 


3.2 - 5.0 


4.1 


26.8 


±2.4 


±4.2 


5.0 - 8.0 


6.4 


17.2 


± 1.7 


±3.0 


8.0 - 12.7 


10.1 


11.5 


± 1.3 


±2.0 


12.7 - 20.1 


16.0 


6.4 


± 1.0 


± 1.1 


20.1 - 31.8 


25.0 


2.20 


±0.55 


±0.38 


31.8 - 80.0 


50.0 


0.57 


±0.20 


±0.10 



Table 5: Q 2 dependence of the elastic cross section a(^i*p — > J/ijjp). 
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Q 2 




da(ep^eJ/?pX)/d Q 2 


[pb/GeV 2 ] 




[GeV 2 l 

L J 




inclusive 






M x > 10 GeV 


2.8 


269 


± 16 (stat.) 


± 43 (syst.) 


95.5 


± 8.8 (stat.) 


± 15.3 (syst.) 


5.0 


127 


± 8 


±20 


44.4 


±4.3 


±7.1 


9.1 


39.5 


±3.2 


±6.3 


16.4 


±2.0 


±2.6 


17.1 


12.0 


± 1.3 


± 1.9 


4.5 


±0.7 


±0.7 


31.1 


2.28 


±0.40 


±0.36 


1.69 


±0.37 


±0.27 


54.5 


0.35 


±0.12 


±0.06 


0.25 


±0.11 


±0.04 


2 

pi ,i, 




da{ep^eJ/^X)/dp 2 ^ 


[pb/GeV 2 ] 




[GeV 2 ] 

L J 




inclusive 






M x > 10 GeV 


1.4 


177 


± 20 (stat.) 


± 44 (syst.) 


73.7 


±11.6 (stat.) 


± 18.4 (syst.) 


2.6 


180 


± 14 


±29 


51.4 


±7.0 


± 8.2 


4.8 


114 


± 8 


± 18 


30.5 


±3.8 


±4.9 


8.6 


43.4 


±3.5 


±6.9 


16.7 


±2.1 


±2.7 


15.5 


11.5 


± 1.2 


± 1.8 


5.0 


±0.7 


±0.8 


27.0 


2.14 


±0.35 


±0.34 


1.28 


±0.27 


±0.20 








da(ep — > eJ/i/jX)/d z 


[Pb] 




z 




inclusive 






M x > 10 GeV 


0.275 


690 


± 240 (stat.) 


± 110 (syst.) 


660 


± 230 (stat.) 


± 110 (syst.) 


0.425 


460 


± 110 


±70 


420 


± 130 


±70 


0.575 


590 


± 80 


± 110 


550 


± 80 


± 100 


0.725 


590 


± 80 


± 160 


470 


± 70 


± 130 


0.850 


950 


± 150 


±310 


730 


± 130 


±240 


0.950 


8350 


±390 


±2000 


1170 


± 110 


±280 


y 




da(ep^eJ/ipX)/dy* 


[Pb] 






inclusive 






M x > 10 GeV 


1.77 


73 


± 25 (stat.) 


± 17 (syst.) 


71 


± 24 (stat.) 


± 16 (syst.) 


2.24 


396 


±38 


±91 


277 


± 31 


± 64 


2.68 


737 


±53 


± 133 


343 


± 35 


± 62 


3.00 


746 


±51 


± 134 


244 


± 28 


±44 


3.32 


702 


±49 


± 126 


191 


±25 


± 34 


3.74 


568 


±62 


± 102 


115 


± 23 


±21 


W 




da(ep^eJ/ipX)/dW 


[pb/GeV] 




[GeV 

L J 




inclusive 






M x > 10 GeV 


50 


11.1 


±1.1 (stat.) 


± 2.0 (syst.) 


4.3 


±0.8 (stat.) 


±0.8 (syst.) 


70 


10.0 


±0.8 


± 1.8 


3.4 


±0.4 


±0.6 


90 


9.4 


±0.7 


± 1.7 


3.9 


±0.4 


±0.7 


110 


9.1 


±0.7 


± 1.6 


3.3 


±0.4 


±0.6 


130 


8.0 


±0.8 


± 1.4 


3.4 


±0.5 


±0.6 


150 


9.0 


± 1.1 


± 1.6 


4.0 


±0.6 


±0.7 


170 


8.5 


± 1.6 


± 1.5 


3.2 


±0.8 


±0.6 



Table 6: Inclusive and inelastic (M x > 10 GeV) differential cross sections for the process 
ep -> e J/tp X in the kinematic region 2 < Q 2 < 80 GeV 2 , 40 < W < 180 GeV and z > 0.2. 
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da(ep^eJ/^X)/d 


z [pb] 




z 




inclusive 






M x > 10 GeV 


0.275 


110 


± 80 (stat.) 


± 20 (syst.) 


110 


± 80 (stat.) 


± 20 (syst.) 


0.425 


210 


± 80 


±30 


190 


±70 


±30 


0.575 


180 


±50 


±30 


160 


± 50 


±30 


0.725 


220 


±50 


±60 


200 


± 50 


±50 


0.850 


250 


±70 


± 80 


190 


± 70 


±60 


0.950 


3750 


±240 


±900 


530 


± 70 


± 130 






da(ep^ eJ/tpX)/dy* [pb] 






inclusive 






M x > 10 GeV 


1.77 


15 


± 10 (stat.) 


± 3 (syst.) 


15 


± 10 (stat.) 


± 3 (syst.) 


2.24 


155 


±22 


±36 


95 


± 17 


±22 


2.68 


328 


±33 


±59 


126 


± 20 


±23 


3.00 


317 


± 32 


±57 


80 


± 15 


± 14 


3.32 


309 


± 34 


±56 


66 


± 14 


± 12 


3.74 


169 


±35 


±30 


31 


± 13 


±6 


W 




da{ep^eJ/ipX)/dW 


[pb/GeV] 




[GeV 2 ] 




inclusive 






M x > 10 GeV 


50 


3.8 


± 0.6 (stat.) 


± 0.7 (syst.) 


1.32 


± 0.42 (stat.) 


± 0.24 (syst.) 


70 


4.2 


±0.5 


±0.8 


1.12 


±0.25 


±0.20 


90 


3.9 


±0.4 


±0.7 


1.32 


±0.23 


±0.24 


110 


4.1 


±0.5 


±0.7 


1.30 


±0.24 


±0.23 


130 


3.2 


±0.5 


±0.6 


1.12 


±0.25 


±0.20 


150 


3.5 


±0.6 


±0.6 


1.63 


±0.35 


±0.29 


170 


2.1 


±0.6 


±0.4 


0.80 


±0.33 


±0.14 



Table 7: Inclusive and inelastic (M x > 10 GeV) differential cross sections for the process 
ep -> e J/ipX in the kinematic region 4 < Q 2 < 80 GeV 2 , p 2 ^ > 4 GeV 2 , 40 < W < 180 GeV 
and z > 0.2. 



w 






<?{l*P - J MX) 


[nb] 




[GeV] 




inclusive 






z < 0.9 




50 


33.0 


± 3.2 (stat.) 


± 5.9 (syst.) 


16.0 


± 3.2 (stat.) 


± 2.9 (syst.) 


70 


43.1 


±3.3 


±7.8 


12.9 


± 1.9 


±2.3 


90 


53.8 


±4.0 


±9.7 


17.7 


±2.3 


± 3.2 


110 


66.3 


±5.2 


± 11.9 


22.1 


±3.1 


±4.0 


130 


72.8 


±7.1 


± 13.1 


28.0 


±4.5 


±5.0 


150 


101 


± 12 


± 18 


33.6 


±6.0 


±6.0 


170 


115 


± 22 


±21 


27.9 


±7.2 


±5.0 



Table 8: W dependence of the inclusive and the inelastic (z < 0.9) cross sections u(YP 
J/ip X) in the kinematic region 2 < Q 2 < 80 GeV 2 , 40 < W < 180 GeV and z > 0.2. 
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